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ABSTRACT 


(Distribution  Limitation  Statement  A) 


The  electron  collision  frequency  is  a  critical  parameter  in  determining  the 
attenuation  of  electromagnetic  signals  transmitted  through  the  plasma  sheath 
surrounding  advanced  reentry  vehicles.  This  study  improves  current  methods  of 
obtaining  this  parameter  for  multicomponent  gas  mixtures.  A  multicomponent 
collision  frequency  model  is  defined  with  temperature  and  species-dependent 
electron  collision  cross  sections.  A  parametric  study  of  collision  frequency  in 
a  gaseous  mixture  of  phenolic  carbon  and  air  Is  performed  using  this  model,  and 
the  results  are  compared  to  calculations  made  using  a  constant  cross-section, 
clean  air  model.  A  study  of  collision  frequency  and  signal  attenuation  in 
turbulent  boundary  layers  is  performed  using  the  multicomponent  and  constant 
c’'oss-section  models.  This  study  shows  that  the  multicomponent  model  can  predict 
collision  frequencies  in  an  ablation  air  mixture  that  are  signif’cantly  different 
frren  those  pr^icted  by  constant  cross-section  models.  These  differences  can 
result  fran  large  cross-section  ablation  molecules  or  from  the  temperature 
dependence  of  the  multicomponent  cross  sections.  The  effect  these  frequency 
differences  have  upon  the  signal  attenuation  are  dependent  on  the  electron  den¬ 
sity  profiles.  The  attenuation  predicted  by  the  two  models  may  or  may  not  differ 
significantly. 
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SECTION  I 
INTRODUCTION 

Hypersonic  reentry  vehicles  encounter  severe  aerodynamic  heating  rates  that 
necessitate  ablative  heatshield  thermal  protection  to  ensure  vehicle  survival- 
Unfortunately,  the  addition  of  heatshield  material  to  the  flow  field  surrounding 
a  reentry  vehicle  may  create  a  plasma  that  can  seriously  attenuate  signals 
tran.;mitted  to  or  from  the  vehicle.  This  signal  attenuation  at  radio  frequencies 
can  be  principally  attributed  to  reflection  and  absorption  of  the  signal  by  free 
electrons  in  the  flow  field. 

Heatshield  materials  are  contaninated  by  small  amounts  of  easily  ionized 
chemical  species,  usually  predominately  sodium,  that  enter  the  flow  field  through 
ablation.  Typically  at  altitudes  below  100,000  feet,  these  ablation  contaminates 
and  the  free  electrons  they  produce  are  confined  to  the  chemically  reacting 
turbulent  boundary  layer  that  enshrouds  the  hypersonic  vehicle.  The  plasma 
sheaths  are  typically  collision-dominated  since  the  charge  particle  densities  are 
usua.ly  several  orders  of  magnitude  smaller  than  the  neutral  particle  densities. 

At  a  given  signal  frequency,  the  attenuation  of  a  plasma  sheath  is  charac¬ 
terized  by  electron  den‘:ity  and  electron  collision-frequency  profiles.  Sophis¬ 
ticated  computer  codes  can  be  used  to  determine  chemical  species  and  electron 
concentrations  in  chemically  reacting  turbulent  boundary  layers.  Present  methods 
fcr  calculating  electron  collision  frequencies  in  the  plasma  sheath  formed  on 
advanced  hypersonic  reentry  vehicles  are  based  on  collision  cross  sections  for 
high-temperature  air  (refs.  1  and  2).  Since  the  ionized  flow  about  these  vehicles 
contains,  in  addition  to  air  species,  a  number  of  ablation  species  with  widely 
varying,  temperatu»'e-deperident  collision  cross  sections,  the  present  electron 
collision  frequency  model  may  be  inadequate  for  making  accurate  predictions  of 
the  plasma  sheath  electromagnetic  signal  attenuations. 

This  study  defines  an  improved  collison  frequency  model.  The  improved  and 
present  collision  frequency  models  are  compared  in  an  Investigation  of  the  effect 
of  ablation  products  in  an  air  mixture  on  electron-neutral  particle  collision 
frequency.  The  two  models  are  compared  for  a  closed  system  parametric  study  and 
typical  reentry  vehicle  boundary  layers.  Signal  attenuation  is  calculated  using 
the  two  models  ^or  e  typical  reentry  vehicle. 
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The  present  collision  frequency  model  employed  in  the  Air  Fo»ce  Weapons 
Laboratory  (AFWL)  nonsimilar  turbulent  boundary  layer  code,  BLIMP  (ref.  3),  ano 
the  AFWL  similar  turbulent  boundary  layer  code,  REBOUND  (ref.  1),  defines  the 
electron-neutral  collision  cross  section  as  a  constant  value  representative  of 
high- temperature  clean  air.  Collision  frequency  is  calculated  with  this  one 
value  of  collision  cross  section,  whereas  the  improved  model  defines  the  colli¬ 
sion  cross  section  as  a  function  of  temperature  and  the  chemical  species. 
Collision  frequency  is  determined  as  the  total  of  the  collision  frequencies  for 
each  specie.  The  improved  model  is  thus  physically  more  realistic  than  the 
present  model . 


2 


AFWL-TR-72-71 


SECTION  II 
METHODOLOGY 


1.  SIGNAL  A1TENUATION  CALCULATION 

To  a  first  approximation,  the  attenuation  of  an  electromagnetic  signal  on 
transmission  through  a  plasma  sheath  can  be  calculated  using  a  plane  wave 
analysis  (ref.  4).  Plasma  property  variations  are  represented  by  dividing  a 
plasma  sheath  into  many  slabs  such  that  the  properties  of  each  slab  can  be 
assumed  homogeneous.  Each  homogeneous  slab  is  represented  by  its  cwnplex 
dielectric  constant 


where  Wp  *  56,000.yii^.  Thus,  transmission  through  a  slab  is  characterized  by 
electron  density  (Ng),  collision  frequency  (v^),  and  the  signal  frequency 
The  propagation  constant  for  each  slab  is  given  by 


•p 


cos  6, 


cos  e. 


(2) 


For  transverse  waves  the  transmitted  electric  and  magnetic  fields  are  related 
to  the  input  electric  and  magnetic  fields  of  a  single  slab  by  the  matrix 
equation 


■e,- 

cosh  Vpd  -  jZp  sinh  Ypd 

■Et' 

i  sinh  y  d 

- - — —  cosh  Y:d  • 

- 1 

_ _ _ _ 1 

3 


AFHL-TR-7'?-71 


or 


'h' 

■Et' 

=  M 

"t 

(4) 


where 


is  the  characteristic  impedance  of  the  slab  given  by 


“P 


(5) 


The  input  and  output  fields  for  n  homogeneous  slabs  are  related,  to  a  close 
approximation,  by  the  product  of  the  n  transmission  matrices  for  the  slabs 
(ref.  4) 


(6) 


The  total  signal  attenuation  for  a  plasma  sheath  divided  into  n  slabs  can  be 
determined  from  the  above  equation,  ffeltiplying  the  logarithm  of  the  ratio  of 
the  input  power  to  the  output  power  by  ten  gives  the  total  attenuation  in  units 
of  decibels. 

2.  COLLISION  FREQUENCY  CALCULATION 

The  collision  frequency  In  a  nfjlticon’p^nent  ionizsd  gas  for  electron-molecule, 
electron-atOTi,  a»Wor  electron-ion  interactions  is  calculated  (ref.  5)  by 


/SkTA 

'^c  yj  rnn 


ns 


n^Qe- 
j  J 


(7) 


This  eqaaticn  is  obtained  acscming 

a.  Maxwellian  electron  velocity  distribution 

b.  Particles  (molecules,  atoms,  ions)  are  smooth,  rigid,  elastic  spheres 
not  surrounded  by  force  f if Ids 


4 


AFWL-TR-72-71 


c.  Particles  and  electrons  Interact  only  at  collision, 
cl.  Mixture  .''f  gcses  is  at  a  uniform  steady  state. 

e.  The  mass  and  diameter  of  an  electron  are  negligible  compared  to  that 
of  any  molecule,  atom,  or  ion. 

Equation  6  indicates  that  collision  frequency  is  dependent  on  electron 
velocity  (i.e.,  temperature),  and  the  sunmation  of  the  product  of  each  species 
population  multipl.ad  by  its  collision  cross  section.  "Qe^"  is  the  collision 

«l 

cross  section  for  momentum  transfer  between  an  electron  and  specie  j.  It  is 
dependent  on  electron  temperature  and  the  particular  specie  identity. 

Ablation  material  and  air  may  react  to  ^orm  a  gaseous  mixture  consisting  of 
a  great  number  of  chemical  species.  Inclusion  of  each  species  in  equation  7  may 
not  be  possible  because  of  a  lack  of  collision  cross-section  data  or  excessive 
numbers  of  species.  The  product  n-Oe-  in  equation  7  may  be  used  to  estimate  the 
relative  contribution  of  an  individual  specie  to  the  total  collision  frequency. 

A  number  of  species  may  be  neglected  in  the  collision  frequency  calculation,  if 
the  term  the  neglected  species  is  small  compared  to  ^n-Qe-  for 

the  remaining  species. 

The  equation  presently  used  for  calculating  collision  frequency  in  an  abla¬ 
tion  contaminated  boundary  layer  is  (ref.  1) 

_  5.814  X  lO^^P 

V - -  (8) 

All  ti.c  assumptions  made  for  equation  7  apply  to  equation  8.  Equation  8  can  be 
derived  frot.i  equation  7  with  the  additional  assumption  that  the  mixture  is  a 
single  component  ideal  gas  with  a  constant  collision  cross  sect'on  which  is 
representative  of  high- temperature  air.  A  molecular  diameter  of  3.65  x  10~®  cm 
(10.4  X  10"’^  cm^  cross  section)  and  a  molecular  we.ght  of  28.962  are  used  for 
this  model. 

3.  ELECTRON-PARTICLE  COLLISION  CROSS  SECTION 

The  probability  that  an  elastic  interaction  will  occur  between  an  electron 
and  a  particle  (atom,  ion,  or  molecule)  is  defined  (ref.  5)  as  the  total  elastic 
scattering  collision  cross  section,  Q^.  The  probability  that  an  electron  will 
be  scattered  into  a  differential  solid  angle,  dn  =  sinededi^,  in  an  elastic 
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collision  is  p(e)  sineded^.  Defining  1(e)  =  p(6)Q^,  gives  the  differential 
scattering  cross  section,  1(6)  sineded<),  and  the  total  cross  section  for 
scattering. 


Qc 


1(6)  sin  ded4» 


(9) 


The  energy  lost  in  a  niomentum  transfer  collision  by  an  electron  scattered 
through  an  angle  6  is  2(l-cos6)"’®/mj.  Accordingly  the  collision  cross  section 
for  momentum  transfer  is  defined  as 


IT  Z-n 

1(6)  (l-cose)  sinededi  (10) 

0  0 

If  1(6)  is  independent  of  e,  then  the  right  hand  sides  of  equations  9  and 
can  be  integrated  to  show  that  =  0^^^.  This  equality  is  approximately  true 
for  slew  electrons  (ref.  7).  The  electron  velocities  encountered  by  a  typical 
reentry  ‘-.icle  are  sufficiently  slow  to  assume  =  0^  because  thermal  electron 
energies  of  less  than  1  electron  volt  correspond  to  a  typical  reentry  vehicle 
temperature  range  of  1000°R  to  lOjCOO^R.  With  this  justification  elastic 
scattering  collision  cross-section  data  were  substituted  for  momentum  transfer 
cross-section  data  when  the  latter  were  unavailable. 

4.  LITERATURE  SURVEY 

A  literature  survey  was  performed  to  collect  cross-section  data  for  chemical 
species  which  might  be  found  in  an  ablation  mixture  of  heatshielo  material  and 
e’r.  Since  phenolic  carbon  and  phenolic  silica  are  typical  heatshield  materials, 
docifflients  containing  cross-sectico  data  fo**  molecules  containing  any  combination 
of  the  elements  C,  H,  N,  0,  A,  and  Si  were  collected  (rets.  8  to  24).  These 
cross-section  data  are  also  applicable  to  the  AFWL  rocket  sled  plasma  experiment. 
The  nose  cone  heatshield  on  the  rocket  sled  is  avcoat  8024  (epoxy  urethane, 

C567  ^950  ®113  ^35^  seeded  with  20  percent  by  weight  cesium  nitrate.  Collision 
data  for  the  heatshield  contaminants,  sodium  and  cesium,  were  also  collected. 

Severs!  sources  of  air  collision  data  were  located,  but  cross-section  data 
for  molecules  containing  silica,  hydrogen,  and/or  carbon  was  limited.  Shkarofsky, 
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Backynski,  and  Johnston  (ref.  8)  have  collected  and  evaluated  data  from  several 
sources  for  each  of  the  air  molecules  N,  NO,  N2,  0,  O2,  and  A.  This  report 
guided  the  selection  of  air  collision  cross-section  data  used  in  the  present 
study.  Cross-section  data  for  the  following  molecules  containing  hydrogen  cr 
carbon  were  collected:  H,  H2,  H2O,  HCN,  C,  CO,  CO2,  C2Hj,,  C2H2  (refs.  6,  9,  18, 
19,  20,  23).  Cesium  and  sodium  crofs-section  data  were  available  in  reference 
6  and  9.  Onl’  el'^mental  silica  cross-section  data  were  found.  The  Bibliography 
of  Low-Energy  Electron  Collision  Cross-Section  Data  (refs.  10  and  11),  a 
continuously  updated  guide  to  material  on  electron  collisions,  contained  the  only 
references  found  for  the  silica  atom.  Because  data  for  molecules  containing 
silica  were  not  found,  phenolic  silica  was  not  investigated  in  this  study.  The 
collected  collision  data  are  presented  in  figure  1. 

Two  sources  were  found  that  contain  a  considerable  amount  of  collision  data. 
Massey  and  Burhop  (ref.  9)  have  compiled  perhaps  the  most  extensive  set  of 
collision  cross-section  data  in  their  text.  Electronic  and  Sonic  Impact  Phenomena. 
Brown  (ref.  6)  used  a  computer  information  retrieval  system  to  compile  cross- 
section  data  in  Basic  Data  of  Plasma  Physics,  1966. 
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SECTION  III 
PARAMETRIC  STUDY 

A  closed  system  parametric  study  was  performed  for  a  uniform  gaseous  mixture 
of  phenolic  carbon  (C^HgO)  and  air  in  chemical  equilibrium.  Temperature,  pres¬ 
sure,  and  mass  fraction  of  ablation  material  in  the  mixture  were  selected  as  the 
paraireters  of  interest  for  this  study.  Parameter  ranges  were  chosen  to  give 
system  states  typical  of  reentry  boundary  layers.  System  pressure  was  held  con¬ 
stant  at  values  of  0.1,  0.5,  1.0,  3.0,  and  5.0  atmospheres,  while  collision 
frequency  calculations  were  made  varying  temperature  and  the  mass  fraction  of 
ablator.  The  temperature  was  varied  from  1045“R  to  10,450“R  (corresponding  to 
electron  energies  of  0.05  eV  to  0.50  eV)  and  the  mass  fraction  of  ablator  was 
varied  from  0  to  50  percent. 

The  selection  of  chemical  species  to  model  the  mixture  was  based  on  availa- 
bilty  of  cross-section  data  and  particle  concentrations  in  typical  reentry 
boundary  layers  predicted  by  computer  using  a  60  species  chemistry  model.  With 
reference  to  the  discussion  in  section  2.2,  the  terms  njQej  were  used  to  determine 
which  chemical  species  could  be  neglected  in  the  collision  frequency  calculations. 
Since  cross-section  data  were  not  available  for  every  species,  an  upper  limit  for 
collision  cross  section  was  assumed  equal  to  the  largest  cross  section  found  in 
the  literature.  The  largest  cross  section  found  was  6020  x  10*^®  cm^  for  the 
large  dipolar  molecule  CH3NO2  at  295®K  (ref.  18).  This  upper  limit  is  approxi¬ 
mately  three  orders  of  magnitude  larger  than  the  cross  section  of  the  density 
dominant  Nj  molecule  in  a  turbulent  boundary  layer.  To  limit  the  magnitude  of 
the  term  rejected  species  to  a  small  fraction  of  the  term 

X,  n.Qe.  for  the  considered  species,  only  species  having  particle  densities 

w  %2 

five  or  more  orders  of  magnitude  smaller  than  the  density  of  N2  were  eliminated 
from  the  collision  frequency  calculations. 

Collision  cross  sections  for  species  considered  in  this  study  are  presented 
in  figure  1.  Air  species  cross  sections  range  from  1.7  x  10"^*  cm^  to  1.4  x 
10*^^  cm^,  approximately  an  order  of  magnitude  of  variation  with  the  exception 
of  the  small  species  argon.  The  species  HCN,  H2O,  CO2,  and  H  have  large  cross 
sections  compared  to  air  species,  therefore  small  concentrations  of  these  mole¬ 
cules  in  an  air  mixture  can  significantly  affect  the  collision  frequency  of  the 
mixture. 
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Collision  frequencies  calculated  for  this  parametric  study  are  presented  in 
figures  2  to  6.  Each  figure  presents  calculations  at  a  single  constant  pressure, 
and  each  curve  depicts  collision  frequency  at  a  particular  constant  ablator  mass 
fraction  as  a  function  of  temperature.  These  curves  appear  similar  in  shape 
from  one  pressure  to  the  next.  The  curves  shift  upward  in  collision  frequency 
as  pressure  is  increased.  Perhaps  the  most  startling  result  of  this  study  is 
the  erratic  variation  of  collision  frequency  with  the  mass  fraction  of  phenolic 
carbon  in  the  contaminated  air  mixture.  This  variation  is  further  illustrated 
in  figure  7  where  collision  frequency  is  plotted  at  constant  pressure  and  tempera 
ture  to  give  a  mass  fraction  dependent  graph. 

Particle  concentrations  in  the  gaseous  mixture  and  corresponding  electron 
collision  cross  sections  can  explain  the  collision  frequency  variations  in 
figures  2  to  7.  Particle  concentrations  are  plotted  for  mixtures  at  pressures 
of  1  atmosphere  and  5  atmospheres  for  temperatures  of  3135°R  and  8360®R  in 
figures  8  to  11.  The  collision  cross  sections  of  figure  1  are  used  in  conjunc¬ 
tion  with  these  particle  concentrations  to  illustrate  the  dependence  of  collision 
frequency  on  these  variables. 

The  curves  of  particle  density  versus  ablator  mass  fraction  (figures  8  to  11/ 
shift  upward  with  increased  pressure,  but  appear  similar  in  shape  for  pressures 
of  1  and  5  atmospheres  at  constant  temperature.  This  upward  shift  and  similarity 
plus  the  assumption  that  collision  cross  sections  are  not  pressure  dependent 
explain  the  upward  shift  and  the  similarity  of  the  collision  frequency  curves. 
However,  the  shape  of  the  curves  in  figures  2  to  6  are  not  identical.  The  fact 
that  the  particle  mole  fractions  are  weakly  pressure  dependent  can  be  seen  in 
figures  8  to  11  that  shew  slight  variations  in  curve  shapes  at  constant  tempera¬ 
ture  with  variation  in  pressure.  These  slight  variations  are  reflected  in  the 
collision-frequency  figures.  For  instance,  the  frequency  curves  for  mass  frac¬ 
tions  of  40  and  50  p-ercent  have  maxima  with  peaks  that  shift  from  approximately 
6500“R  to  7800°R  as  pressure  is  increased  from  0.1  atmospheres  to  5  atmospheres. 
Figure  12  is  presented  to  show  that  this  shift  is  due  to  a  shift  in  the  mole 
fraction  of  the  large  cross-section  molecule  HCN. 


The  erratic  variation  of  collision  frequency  with  ablator  mass  fraction  can 
be  related  to  the  particle  concentrations  of  three  large  cross-sections  molecules 
HCN,  HjO,  and  CO2.  Figiire  7  shows  two  maxima  in  the  collision  frequency  curves 
at  a  temperature  of  3135®?..  Figures  8  and  10  show  that  the  lesser  of  the  two 
peaks  corresponds  to  peaks  in  tne  concentrations  of  HjO  and  COj,  while  the  larger 
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Figure  3.  Parametric  Study — Collision  Frequency  at  0.5  Atirospheres  Pressure 
with  Mass  Fraction  of  Ablator  Varying  from  0  to  50  Percent 
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Figure  4.  Parametric  Study — Collision  Frequency  ai  1.0  Atmospheres  Pressure 
with  Mass  Fraction  of  Ablator  Varying  from  0  to  50  Percent 
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Parametric  Study-Collision  Frequency  at  3.0  Atmospheres  Pressure 
with  Mass  Fraction  of  Aolator  Varying  from  0  to  50  Percent 
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Figure  Study- Collision  Frequency  at  5.0  Atmospheres  Pressure 

<ifith  Mass  Fraction  of  Ablator  Varying  from  0  to  50  Percent 
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Figure  7.  Parametric  Study— Collision  Frequency  at  Varying  Ablator 
Mass  Fraction  with  Temperature  and  Pressure  Held  Constant 
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Figure  12.  Parametric  Study— HCN  Concentration  in  a  Mixture  with  an  Ablator 
Mass  Fraction  of  50  Percent  at  Pressures  of  1.0  and  5.0  Atmospheres 
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peak  follows  the  concentration  of  HCN.  At  a  temperature  of  8360°R,  the  density 
of  these  molecules  relative  to  the  total  density  is  smaller  and  the  cross  sections 
are  much  smaller.  Thus,  the  effect  of  HCN,  HjO,  and  CO2  on  collision  frequency 
is  small  and  collision  frequency  does  not  vary  drastically  with  mass  fraction  at 
8360®R.  There  is  a  slight  increase  in  collision  frequency  at  5  atm  and  8360“R. 
This  increase  corresponds  to  a  particle  density  increase  of  HCN.  The  curves  for 
8360°R  are  mainly  the  result  of  superposition  of  collision  frequencies  for  the 
N2,  H,  and  HCN  molecules. 

A  general  trend  of  the  collision  frequency  is  to  decrease  with  increasing 
temperature.  Three  factors  cause  this  decrease.  The  total  number  density 
decreases  with  increasing  temperature  and  the  relative  concentration  of  large 
cross-section  molecules  decreases  with  increasing  temperature.  The  cross  section 
of  the  large  molecules  HCN,  H2O,  CO2,  and  H  decrease  with  increasing  temperature. 
Two  factors  opposing  this  trend  are  the  electron  velocity  which  increases  as  the 
square  root  of  temperature  and  the  slight  increase  of  the  air  species'  cross- 
sections  with  increasing  temperature. 

Collision-frequency  curves  predicted  by  the  constant  cross-section  model  and 
by  the  multicomponent  model  are  compared  in  figures  2  to  7.  By  definition  the 
constant  cross-section  model  cannot  account  for  frequency  variations  caused  by 
specie  concentration  variations  at  a  given  pressure  and  temperature.  Therefore, 
the  constant  cross-section  model  predicts  single  curves  in  figures  2  to  6  and 
straight  lines  in  figure  7  for  varying  ablator  mass  fraction.  These  figures 
show  a  factor  of  five  maximum  deviations  in  the  collision  frequencies  predicted 
by  the  two  models.  Even  for  clean  air  the  two  models  predict  significantly 
different  collision  frequencies.  The  constant  cross-section  model  predicts 
larger  collision  frequencies  because  the  constant  cross  section  is  larger  than 
the  multicomponent  collision  cross  sections.  As  temperature  is  increased  the 
constant  cross  section  and  multicomponent  air  cross  sections  are  more  nearly 
equal . 
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SECTION  IV 

STUDY  OF  TURBULENT  BOUNDARY  LAYERS 

1.  HIGH-VELOCITY  REENTRY  TRAJECTORY;  CASES  1  TO  6 

Collision-frequency  and  signal -attenuation  analyses  were  perforiped  for  six 
points  on  a  typical  hypersonic  reentry  trajectory.  The  reentry  vehicle  studied 
was  a  sphere-cone  with  a  1/2-inch  nose  radius  and  an  8-degree  half-angle  after¬ 
body.  The  afterbody  heatshield  was  phenolic  carbon  with  an  assumed  sodium 
contamination  level  of  100  ppm.  A  nose  cap  material  was  not  assumed  since  the 
solution  method  could  take  into  account  only  one  wall  material.  The  body  station 
studied  was  56.35  streamwise  inches  from  the  stagnation  point.  The  six  trajec¬ 
tory  cases  considered  are  listed  in  table  1. 


Table  I 

REENTRY  TRAJECTORY 


Case 

Altitude 

(feet) 

Velocity 
(ft/ sec) 

1 

74,820 

22,340 

2 

52,530 

21 ,800 

3 

25,930 

20,040 

4 

16,080 

18,870 

5 

11,390 

18,190 

6 

6,870 

17,470 

The  boundary  layer  chemical  species  and  electron  concentration  profiles 
required  for  the  collision  frequency  and  signal  attenuation  analyses  were  obtained 
from  the  AFWL  similar  turbulent  boundary  layer  computer  code  REBOUND,  {ref.  1). 
Boundary  layer  profiles  of  velocity,  enthalpy,  and  ablator  mass  fraction  are 
obtained  by  this  program  using  the  Crocco  similarity  assumptions.  The  boundary 
layer  enthalpy  is  determined  from  a  real  gas  wall  and  edge  enthalpy  applying  the 
Crocco  enthalpy  profile.  Species  concentrations  are  determined  from  pressure, 
enthalpy,  and  ablation  mass  fraction  by  assuming  that  the  gaseous  mixture  is  in 


chemical  equilibrium.  For  a  given  thermodynamic  state*  the  chemical  equilibrium 
composition  is  determined  by  minimization  of  the  Gibb's  free  energy. 

The  results  of  the  boundary  layer  analyses  are  presented  in  figures  13  to  24. 
Collision  frequencies  calculated  using  the  multicomponent  and  constant  cross- 
sectiur  models  are  compared  in  figure  !3  to  18.  Temperature  profiles  are  pre¬ 
sented  in  figure  19  and  ablator  mass  fraction  profiles  are  presented  in  Figure 
20.  Species  concentration  profiles  for  case  3  are  shown  in  figure  21.  The 
electron  density  profiles  for  the  six  cases  are  shown  in  figure  22.  The  boundary 
layer  one-way  signal  attenuation  for  the  six  cases  is  shown  in  figures  23  and  24. 

Both  the  multicomponent  and  constant  cross  section  models  predict  collision 
frequency  profiles  which  appear  similar  for  each  of  the  six  trajectory  cases. 

As  a  general  trend,  the  collision  frequency  curves  shift  upward  with  decreasing 
altitude.  This  is  the  result  of  increased  boundary  layer  density  with  decreased 
altitude.  An  examination  of  the  temperature  and  specie  concentration  figures 
can  explain  the  shape  and  similarity  of  the  curves. 

Recalling  equation  8  and  the  fact  that  pressure  is  constant  across  a  boundary 
layer,  it  is  observed  that  the  constant  cross  section  collision  frequency  is 
inversely  proportional  to  the  square  root  of  the  temperature.  This  relationship 
is  depicted  by  the  collision  frequency  curves  (figures  13  to  18)  and  the  tempera¬ 
ture  curves  (figure  19).  The  frequency  curves  decrease  from  wall  to  approximately 
10"^*  feet  into  the  boundary  layer.  From  10“**  feet  to  the  boundary  layer  edge, 
the  collision  frequency  increases.  Figure  19  shows  that  the  temperature  increases 
frcHu  the  wall  to  approximately  10“**  feet  and  then  decreases  from  10“**  feet  to  the 
boundary-layer  edge. 

Near  the  ablating  wall  the  multicomponent  model  predicts  higher  collision 
frequencies  than  the  constant  cross-section  model,  but  over  a  large  portion  of 
the  boundary  layer  the  multicomponent  model  predicts  lower  collision  frequencies 
than  the  constant  cross-section  model.  This  deviation  is  clarified  by  the  cross 
sections  of  the  species  (figure  1),  the  mass  fraction  profiles  (figure  20),  and 
the  concentration  profiles  for  case  3  (figure  21),  which  are  typical  of  the  six 
trajectory  cases.  The  shapes  of  the  multicomponent  collision  frequency  profiles 
are  dependent  on  temperature  and  species  concentration,  the  pressure  being 
constant.  The  large  cross  section  of  HCN  and  its  number  density  profile,  indicate 
that  this  molecule  makes  a  large  contribution  to  the  collision  frequency  near  the 
wall.  In  conjunction  with  other  significant  contributors  such  as  N2,  H,  and  CO, 
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Figure  15.  Reentry  Trajectory  Case  3 — Boundary  Layer  Collision  Frequency  Profile 


Figure  16.  Reentry  Trajectory  Case  ^--Boundary  Layer  Collision  Frequency  Profile 
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Figure  17.  Reentry  Trajectory  Caoe  5--Boundary  Layer  Collision  Frequency  Profile 
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Figure  18,  Reentry  Trajectory  Case  6--Boundary  Layer  Collision  Frequency  Profile 
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Figure  21.  Reentry  Trajectory  Case  3--Boundary  Layer  Species  Concentration  Profiles 
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Figure  22.  Reentry  Trajectory  Study— Boundary  Layer  Electron  Concentration  Profiles 
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Figure  24.  Reentry  Trajectory  Study — Boundary  Layer  Signal  Attenuation 


AFWL-TR-72-71 


the  HCN  molecule  makes  the  multicomponent  collision  frequencies  greater  than 
the  collision  frequencies  predicted  by  the  constant  cross-section  model  near  the 
wall.  Over  a  largi  portion  of  the  boundary  layer  the  iiwlti component  collison 
frequency  curve  is  relatively  flat,  in  spite  of  an  increase  in  total  particle 
density  towards  the  boundary  layer  edge.  The  mass  fraction  of  ablation  material 
decreases  toward  the  boundary  layer  edge  and  the  large  ablation  species  become 
less  significant.  The  multicomponent  collision  frequencies  become  lower  than 
the  constant  cross-section  frequencies  at  approximately  10"^  feet  because  the 
temperature  decreases  toward  the  boundary  layer  edge  and  causes  a  decrease  in 
the  cross  sections  of  the  clean  air  molecules.  The  slight  rise  in  the  multi- 
component  frequency  curve  near  10“^  feet  corresponds  to  the  densities  of  the 
large  CO2  and  H2O  molecules. 

Electron  densities  that  were  obtained  from  chemical  equilibrium  calculations 
are  presented  in  figure  22.  Although  several  ionized  species  were  included  in 
this  study,  the  free  electrons  predominately  resulted  from  ionized  sodium. 

Because  the  ablative  material  contains  100  ppm  of  sodium  by  weight,  the  mass 
fraction  of  sodium  is  given  by  the  ablator  mass  fraction  profiles  multiplied  by 
io-\ 

Because  the  amount  of  sodium  in  the  boundary  layer  decreases  from  the  wall 
to  the  edge  in  conjunction  with  the  ablator  mass  fraction  {figure  20),  the 
electron  density  might  be  expected  to  decrease  also.  The  electron  density  does 
decrease  toward  the  edge,  but  not  monotonically.  The  electron  density  is  depen¬ 
dent  on  mixture  composition,  pressure,  and  temperature.  It  is  noted  that  both 
peak  temperatures  and  peak  electron  densities  occur  in  the  boundary  layer  region 
between  10’'*  and  10"^  feet. 

Signal  attenuation  calculations  were  made  using  the  collision  frequencies 
plotted  in  figures  13  to  18  and  the  electron  densities  plotted  in  figure  22. 
Comparing  the  multicomponent  signal  attenuation  to  the  constant  cross-section 
signal  attenuation,  a  maximum  deviation  cf  approximately  10  percent  is  observed 
(figures  23  and  24).  This  signal  attenuation  deviation  is  considerably  less 
drastic  than  the  deviations  in  collision  frequency  of  figu**es  15  to  18- 

Equation  1  shows  that  the  complex  dielectric  constant  decreases  as  electron 
density  increases  and  as  collision  frequency  decreases.  The  propagation  constant 
for  the  transmitted  signal  is  proportional  to  the  square  root  of  the  complex 
dielectric  constant.  Thus,  if  the  plasma  sheath  is  divided  into  equal  slabs, 
the  greatest  signal  attenuation  will  occur  in  the  slabs  with  the  largest  ratio 
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of  -jjj — /^iu^  +  v^j.  Observation  of  the  collision  frequency  and  the  electron 

density  curves  indicate  why  the  attenuation  calculations  varied  moderately  for 
the  two  models.  The  largest  deviation  of  collision  frequencies  occurred  in  two 
regions.  The  first  region,  from  the  wall  to  approximately  10"^  feet,  is  very 
thin  rind  thus  had  little  effect  on  the  total  signal  attenuation.  The  second 
region,  from  1.5  x  10"^  feet  to  the  edge,  is  relatively  large  but  the  electron 
density  is  very  small  and  thus  signal  attenuation  is  small.  Most  of  the  signal 
attenuation  occurred  in  the  remaining  region,  from  10"^  feet  to  1.5  x  10"^  feet. 
Figure  25,  which  is  a  linear  plot  of  the  collision  frequencies  and  electron 
density  for  case  3,  clearly  shows  that  most  of  the  signal  attenuation  occurs  in 
this  region.  The  collision  frequency  deviations  in  this  region  are  moderate  and 
explain  the  moderate,  i.e.,  less  than  10  percent,  signal  attenuation  deviations. 

2.  REENTRY  CASE  ANALYZED  BY  BLIMP  CODE;  CASE  ^ 


Because  accurate  determination  of  chemical  species  concentrations  is  impor¬ 
tant  r.o  the  multicomponent  collision  frequency  calculations,  it  is  appropriate 
to  include  a  reentry  case  in  which  the  chemistry  is  completely  coupled  to  the 
boundary  layer  equations.  The  AFWL  nonsimilar  chemically  reacting  turbulent 
boundary  layer  computer  code,  BLIMP  (ref.  3)  solves  these  coupled  equations. 

This  code  assumes  chemical  equilibrium  and  uses  real  gas  properties  in  the  flow 
equations.  The  vehicle  studied  for  this  case  was  a  sphere-cone  with  a  0.238-inch 
nose  radius  and  an  8-degree  half-angle  afterbody.  The  nose  was  graphite,  and  the 
afterbody  heatshield  was  phenolic  carbon.  The  heat-shield  material  was  assumed 
to  contain  50  ppm  of  sodium  by  weight.  A  streamwise  station  located  at  3,5  feet 
from  the  stagnation  point  was  analyzed  at  85,000  feet  altitude  and  22,950  feet 
per  second  velocity. 

Collision  frequency,  specie  concentration,  temperature,  and  electron  density 
profiles  are  presented  for  this  case  in  figures  26  to  29,  respectively.  The 
multicomponent  model  is  seen  to  predict  much  higher  collision  frequencies  than 
the  constant  cross-section  model  in  the  vicinity  of  the  wall,  i.e.,  from  the  wall 
to  5  X  10“**  feet.  In  the  region  from  5  x  10“**  feet  to  1.2  x  10“^  feet  to  the 
edge,  the  collision-frequency  curves  diverge.  An  explanation  of  the  deviation 
between  the  two  models  is  given  by  the  specie  concentration  profiles  in  figure 
27.  The  shape  of  the  constant  cross-section  profile  is  determined  completely  by 
temperature  in  figure  27,  because  pressure  is  constant.  Significant  concen¬ 
trations  of  the  large  HCN,  H,  and  CO  molecules  make  the  specie  dependent 
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Figure  25.  Reentry  Trajectory  Case  3>-Boundary  Layer  Electron 
Concentration  and  Collision  Frequency  Profiles 
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multicomponent  frequency  much  greater  than  the  constant  cross-section  frequency 
in  the  vicinity  of  the  wall.  In  the  region  from  5  x  10”**  feet  to  approximately 
1.2  X  10” feet,  the  electron  collisions  are  predominately  with  H2O,  CO2,  and  Nj. 
From  1.2  x  10“"  feet  to  the  edge,  the  clean  air  N2  and  O2  collisions  are  predomi¬ 
nant-  The  multicomponent  ccl lision-frequency  curve  is  observed  to  correspond  to 
the  concentration  fluctuations  and  temperature  dependent  cross  sections  of  these 
species  in  the  latter  two  regions.  The  divergence  of  the  frequency  curves  near 
the  edge  occurs  because  there  the  constant  cross-section  frequency  is  inversely 
proportional  to  temperature  and  thus  increases  as  the  temperature  decreases 
toward  the  edge.  In  contrast,  the  multicomponent  frequency  decreases  because  the 
concentrations  of  the  large  molecules  decrease  toward  the  edge  and  the  cross 
sections  of  the  clean  air  molecules  decrease  as  the  temperature  decreases  toward 
the  edge. 

Although  the  collision  frequency  models  predict  drastically  different  fre¬ 
quencies  in  the  vicinity  of  the  wall,  the  attenuation  calculations  for  a  2.5  GHz 
transmitted  signal  vary  only  moderately.  The  multicomponent  model  predicted  an 
attenuation  of  5.82  dB,  and  the  constant  cross-section  model  predicted  6.65  dB. 
The  constant  cross-section  value  is  14.3  percent  higher  than  the  imjlticanponent 
value.  An  explanation  for  the  moderate  difference  in  attenuation  calculations 
is  given  by  the  collision  frequency  curves  and  the  electron-density  curve.  As 
in  previous  cases,  the  region  near  the  wall  where  the  collision  frequencies  vary 
drastically  is  very  thin.  The  electron  density  in  this  region  is  not  extranely 
large  so  that  attenuation  in  this  region  is  estiiiated  to  be  small  compared  to 
the  attenuation  in  the  remainder  of  the  boundary  layer.  In  the  region  from 
approximately  10”^  feet  to  the  edge,  the  electron  density  is  very  small  and  the 
attenuation  is  again  small.  Thus  most  of  the  signal  attenuation  occurs  in  the 
region  betvreen  5  x  10”**  feet  and  10~^  feet.  The  collision  frequencies  in  this 
**egior.  are  moderately  different  and  thus  moderately  different  signal  atfenuations 
are  predicted. 

3.  ROCKET  SLED;  CASE  8 

A  final  case  studied  is  the  boundary  layer  formed  about  the  nose  cone  of  a 
hypersonic  rocket  sled.  A  rocket  sled  experiment  was  designed  by  AFWL  to  obtain 
experimental  data  to  help  validate  the  AFWL  nonsimilar  chemically  reacting  turbu¬ 
lent  boundary  layer  code  BLIMP.  The  experiment  also  established  the  feasibility 
of  using  a  rocket  sled  to  form  a  plasma  sheath.  The  sled  nose  cone  is  a  sphere 
cone  with  a  1/2-inch  nose  radius  and  a  9-degrpe  half-angle  cone.  It  has  a 
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tungsten  nos?  tip*  and  the  afterbody  Is  covered  v/ith  Avcoat  3024  (epoxy  urethane, 
^567  ^950  ^113  ^35^  which  is  a  low-perfcrsiance  heatshield  material.  This  material 
Is  seeded  with  20  percent  by  weight  of  CsN’Og,  which  adds  easily  Ionized  cesium 
atoms  to  the  boundary  layer. 

Results  for  this  nonsimilar  boundary  layer  analysis  for  a  streamwise  station 
of  1.696  feet  are  shown  in  figures  30  to  33.  Figure  30  shows  that  the  multi- 
conponent  collisicn  frequency  was  predicted  to  be  much  higher  than  the  constant 
cross-section  collision  frequency  over  a  iarge  portion  of  the  boundary  layer. 

As  in  the  previofis  cases  the  constant  cross-secticn  frequency  profile  cm  be 
directly  related  to  the  temperatu»*e  profile  (figure  31).  Particle  densities 
(figure  32)  in  conjunction  with  cross-section  data  show  that  the  large  CO-, 
and  Cs  molecules  together  with  the  molecule  are  chiefly  responsible  for  the 
behavior  of  the  multicomoonent  collision  frequency.  Figures  30  and  33  show  that 
the  greatest  concentration  of  electrons  in  the  boundary  layer  occurs  m  a  region 
where  the  collision  frequencies  vary  by  a  factor  of  approximately  two.  This 
large  difference  is  reflected  in  the  signal  attenuation  calculation.  For  a  two 
GHz  transmitted  signal,  the  multicomponent  model  predicts  an  attenuation  of  6.80 


dB,  and  the  constant  cross  section  model  predicts  an  attenuation  of  10.6  dB.  The 
constant  cross-section  prediction  is  56  percent  higher  than  the  multicomponent 
prediction. 
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Figure  30.  Rocket  Sled  Analysis--8oundary  Layer  Collision  Frequency  Profile 


Y- NORMAL  DISTAKCE  FROM  WALL  (10^  f« 

Figure  31.  Rocket  Sled  Analysis— Boundary  Layer  Temperature  Profile 


ELECTRON  CONCENTRATION  (IO'<  PARTICLES/CENTIMETER  3) 


AFWL-TR-72-71 


Y- NORMAL  DISTANCE  FROM  WALL  (10'*^  feet) 


Figure  35.  Rocket  Sled  Analysis-Boundary  Layer 
Electron  Concentration  Profile 
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SECTION  V 

CONCLUSIONS  AND  RECOWENOATIONS 

A  parametric  study  showed  that  a  gaseous  mixture  of  phenolic  carbon  and  air 
can  contain  a  large  variety  of  chemical  species  with  widely  varying  collision 
cross  sections.  The  density  and  the  temperature-dependent  cross  sec ton  of  each 
specie  are  taken  into  account  by  the  multicomponent  collision  frequency  model. 
This  model  showed  that  the  collision  frequency  in  an  ablation  air  mixture  can 
vary  radically  as  the  mass  fraction  of  ablator  is  varied.  Since  the  constant 
cross-section  model  defines  collision  frequency  as  a  function  of  only  pressure 
and  temperature,  the  chemical  composition  of  an  ablation  air  mixture  has  no 
effect  on  collision  frequency  predictions  made  using  this  model.  The  parametric 
study  showed  that  collision  frequency  predictions  made  by  the  two  models  can  vary 
significantly.  The  validity  of  representing  the  collision  cross  sections  in  an 
ablation  mixture  by  a  single  constant  cross  section  is  questionable  in  view  of 
the  radical  variations  of  species  and  thus  cross  sections  with  variations  of 
ablstor  mass  fraction.  Because  the  multicwnponent  model  is  physically  more 
realistic,  it  appears  to  be  more  reliable  than  the  constant  cross-section  model 
for  analyzing  such  mixtures. 

Significantly  different  collision  frequencies  in  a  plasma  can  be  predicted 
by  the  two  models.  Near  the  vehicle  wall,  large  ablation  molecules  can  make  the 
multicomponent  collision  frequency  much  larger  than  the  constant  cross-section 
collision  frequency.  The  region  of  the  boundary  layer  in  which  the  ablation 
species  are  important  is  dependent  cr.  ablation  rate,  molecular  species  present, 
and  the  thermodynamic  state.  Near  the  boundary  layer  edge  the  ablation  species 
concentrations  are  low  enough  to  be  negligible.  Here  the  constant  cross  section 
and  nwlti component  models  can  predict  significantly  different  collision  fre¬ 
quencies.  The  cross  sections  of  the  clean  air  species  decrease  with  decreasing 
temperature.  Since  the  constant  cross-section  model  predicts  that  collision 
frequency  is  inversely  proportional  to  the  square  root  of  temperature,  the  col¬ 
lision  frequencies  predicted  by  the  two  models  can  diverge  in  a  region,  such  as 
toward  the  boundary-layer  edge,  »d>ere  temperature  decreases. 

Predictions  of  the  transmitted  signal  attenuation  using  the  two  models  can 
vary  significantly.  For  the  six  trajectory  cases  examined,  the  multicomponent 
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model  predicted  higher  attenuations  than  the  constant  cross-section  model .  The 
constant  cross-section  model  predicted  a  collision  frequency  that  was  14.3 
percent  higher  than  the  multicomponent  model  for  the  85,000  feet  altitude  coupled 
chemistry  case.  The  rocket  sled  case  showed  the  largest  difference  in  attenua¬ 
tion  predictions.  Because  the  H2O,  CO2,  and  Cs  molecules  caused  large  collision 
frequency  deviations  in  a  large  region  of  the  boundary  layer  and  because  most  of 
the  free  electrons  were  concentrated  in  this  region,  a  large  deviation  ir  attenua¬ 
tion  predictions  occurred.  The  constant  cross-secti~  .  model  predicted  ar 
attenuation  that  was  56  percent  higher  than  the  multicomponent  attenuation. 

These  cases  have  shown  that  the  multicomponent  attenuation  can  be  significantly 
higher  or  lower  than  the  constant  cross-section  attenuation.  These  cases  have 
also  shown  that  the  variations  in  attenuation  predictions  can  be  less  severe 
than  the  variations  in  collision  frequency  predictions  depending  on  the  plasma 
electron  density  profile. 

The  six  trajectory  cases  have  shown  that  significant  concentrations  of  the 
large  ablation  molecules  may  be  restricted  to  an  extremely  sinall  region  of  a 
boundary  layer.  In  cases  similar  to  these,  multicomponent  collision  frequency 
predictions  based  on  clean  air  cross  sections  only  may  be  sufficient  to  predict 
signal  attenuation.  At  temperatures  from  approximately  8000‘’R  to  10,CCu‘^h  the 
constant  cross  section  which  is  used  for  the  constant  cross-section  model  appears 
to  represent  the  average  cross  section  of  clean  air  with  reasonable  accuracy. 

Use  of  the  constant  cross-section  model  to  predict  collision  frequencies  at  these 
temperatures  should  not  cause  significant  error.  It  may  be  expected  that  the 
contribution  of  large  ablation  molecules  to  the  total  electron  collision  fre¬ 
quency  will  increase  with  increasing  ablation  rate.  Since  the  alkali  metal 
contaminates  in  a  heatshield  material  are  the  main  contributors  of  free  electrons 
in  a  boundary  layer  plasma,  the  electron  density  may  be  expected  to  increase  with 
increasing  ablation  rate.  It  may  thus  be  conjectured  that  as  a  signal  attenua¬ 
tion  increases  because  of  increasing  electron  density,  the  contribution  of  abla¬ 
tion  species  to  the  total  collision  frequency  will  be  increasingly  important, 
and  thus,  it  becomes  increasingly  important  to  use  the  multicomponent  model  to 
obtain  accurate  predictions.  An  inspection  of  the  product  n^  Oej  for  species  in 
a  boundary  layer  can  be  used  to  estimate  how  much  detail  to  include  in  a  collision 
frequency  analysis. 
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Use  of  the  multicomponent  collision  frequency  model  is  restricted  by  two 
major  factors:  1)  collision  cross  sections  must  be  available  for  important 
chemical  species,  and  2)  an  accurate  chemical  composition  analysis  is  required. 

Only  a  phenolic  carbon  ablator  has  been  investigated  in  this  study.  Of  equal 
importance  are  mixtures  of  phenolic  silica  and  air.  Neither  analytical  or  experi¬ 
mental  values  have  been  found  in  the  open  literature  for  many  of  the  ablation 
species  that  contain  silica.  It  is  suspected  that  collision  cross-section  data 
do  not  exist  for  many  of  these  molecules.  Therefore,  fundamental  analyses  or 
experimentation  to  obtain  cross-section  data  should  be  performed  if  this  model 
is  to  be  applied  to  silica  materials.  The  multicomponent  model  is  dependent  on 
species  concentrations  of  a  gaseous  mixture.  It  has  been  shown  that  small 
amounts  of  large  molecules  can  greatly  affect  the  collision  frequency.  Thus,  an 
accurate  analysis  of  the  chemical  composition  of  a  gaseous  mixture  must  be 
obtained  before  this  model  can  predict  a  reliable  and  accurate  collision  frequency. 
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